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Abstract 

We study the thermal-noise spectrum of multi-loop pendulum suspensions for test masses in 
interferometric gravitational-wave detectors. The dependence of the thermal noise on suspension 
parameters and on properties of the wire material is discussed for the situation in which the losses 
are dominated by the internal friction in the pendulum wires. 

PACS: 04.80.Nn; 05.40. Ca; 62.40.+i 

1 Introduction 

The thermal noise is expected to be one of the main limitations on the sensitivity of long-baseline 
interferometric gravitational- wave detectors like LIGO and VIRGO (l|, [2j. Thermal fluctuations of 
internal modes of the interferometer's test masses and of suspension modes will dominate the noise 
spectrum at the important frequency range between 50 and 200 Hz (seismic noise and photon shot 
noise dominate for lower and higher frequencies, respectively). The thermal fluctuations in pendulum 
suspensions were studied both theoretically and experimentally in a number of works (see, e.g. Refs. ^, 
@) i> El @> [| 10, 11, 12, 14, 15]). The predictions of the thermal-noise spectrum in interferometric 



gravitational-wave detectors combine theoretical models (with the fluctuation-dissipation theorem of 



statistical mechanics [16] serving as a basis) and experimental measurements of quality factors of 



systems and materials involved. It is usually assumed that losses in the suspensions will occur mainly 



due to the internal friction in the wires, which is related to anelasticity effects [17, 18]. This assumption 
will be correct only provided that all the losses due to interactions with the external world (friction 
in the residual gas, dumping by eddy currents, recoil losses into the seismic isolation system, friction 
in the suspension clamps, etc.) are made insignificant by careful experimental design. 

In the present work we consider a multi-loop pendulum suspension and study the dependence of 
the thermal-noise spectrum on properties of the wire material and on suspension parameters. The 
thermal-noise spectral density x 2 (u>) depends strongly on the type of the internal friction in the wires. 
We consider two possibilities: (i) the wire internal friction with a constant loss function and (ii) the 
thermoelastic damping mechanism [jOJ. The main conclusion is that the thermal noise can be reduced 
by increasing the number of suspension wires, especially in the case of the thermoelastic damping. 
This conclusion is valid as long as the dissipation due to the friction in the suspension clamps is 
insignificant. 
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2 Thermal-noise spectrum for a pendulum suspension 



In interferometric gravitational-wave detectors, the test masses are suspended as pendulums by one 
or two loops of thin wires. We will consider a multi-loop suspension with the wires attached to the 
bob near the horizontal plane which cuts the bob through its center of mass. We will also assume that 
the mass of the wires is much smaller than the mass of the bob. In such a multi-loop suspension the 
rocking motion of the test mass is essentially suppressed and the main contribution to the thermal- 
noise spectrum is due to the pendulum mode and the violin modes. Then one can write the suspension 
thermal-noise spectral density as a sum, 

x 2 (u;) = xl(iv) + x 2 v (u J ), (1) 

of the pendulum-mode contribution, x 2 (lo), and of the violin-modes contribution, x^(w). 

According to the fluctuation-dissipation theorem, the pendulum-mode contribution can be ex- 
pressed as H 

2 4k B T Wp^pH 



xJu) 



^ ' COM (u;2_ w 2)2 +a; 402' V , 

where ks is Boltzmann's constant, T is the temperature, M is the pendulum mass, </> p (w) is the loss 
function, uj v = (g/L) 1 / 2 is the pendulum frequency, g is the acceleration due to the Earth gravity field, 
and L is the pendulum length. Note that the spectral density x 2 {uj) is written explicitly in terms of 
the angular frequency u, but in fact the density is with respect to the linear frequency / = uj/2tt and 
x 2 {oj) is measured in units of m 2 /Hz. 

The loss function is a measure of the energy dissipation. Let S be the total energy of a dissipative 
oscillator (assuming that the losses are small) and AS be the energy dissipated per cycle. Then 

The energy of the pendulum consists of two parts: the gravitational energy S gr and the elastic energy 
S e i due to the bending of the wire. The gravitational energy is lossless; provided that all the losses due 
to interactions with the external world are made insignificant by careful experimental design, the as- 
sumption is made that the losses are dominated by internal friction in the wire material. Consequently, 
AS = AS e \, and one obtains 

4>p = £ P 0w, (4) 

where </> w = AS c \/(2irS c \) is the loss function for the wire itself which occurs due to anelastic effects 
in the wire material, and £ p = (£ e i/£gr)p is the ratio between the elastic energy and the gravitational 
energy for the pendulum mode. The elastic energy depends on how many wires are used and how they 
are attached to the pendulum bob. In the multi-loop configuration we consider, the wires bend both 
at the top and the bottom, so £ p ~ (fceL)™ 1 , where k^ 1 ~ (EI/T) 1 / 2 is the characteristic distance 
scale over which the bending occurs. Here, T is the tension force in the wire, E is the Young modulus 
of the wire material, and / is the moment of inertia of the wire cross section (/ = ^vrr 4 for a cylindrical 
wire of radius r). For a suspension with N wires (the number of wires is twice the number of loops), 
T = Mg/N, and one obtains 

NVTEI _ 1 I EIN 
MgL "lVW () 
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For LIGO suspensions, f p = uj p /2tt is about 1 Hz. This is much below the working frequency range 
(near 100 Hz), so we may assume uj p /uj <C 1. Also, the loss function is very small, (j) p < 10~ 5 . Then 
the pendulum-mode contribution to the thermal noise spectrum is 



" — mJ> = "P- V M 3 ^ ■ (6) 

The contribution of the violin modes to the thermal noise spectrum is given by 



where n = 1, 2, 3, . . . is the mode number. The angular frequency of the nth mode is 



mr 



2 1 / nvr x 2 



1 + — + 



(8) 



where p is the linear mass density of the wire. For heavily loaded thin wires like in LIGO, k e <S L, 

so 

n-K Ft 

"—TV V (9) 

This is just the angular frequency of the nth transverse vibrational mode of an ideal spring. The 
effective mass of the nth violin mode is 

where we took expression (Joj) for u n and T = Mg/N. This effective mass arises because the violin 
vibrations of the wire cause only a tiny recoil of the test mass M. The loss function for the nth violin 
mode is 

<t>n = £n0w, (11) 

where £ n = (£ e i/£ gr )n is the ratio between the elastic energy and the gravitational energy. This ratio 
is[§ 

Since k e L 3> 1, for first several modes the energy ratio is approximately 



2 EIN , , 

= IV W (13) 

This expression takes into account only the contribution to the elastic energy due to wire bending 
near the top and the bottom. For higher violin modes, one should also consider the contribution due 
to wire bending along its length, which leads to Eq. (|l^). 

Typical values of f\ = uji/2it are from 250 to 500 Hz. If we are interested in the thermal spectral 
density near 100 Hz, we can assume u> 2 <C w 2 . Then we have approximately 

Sk B Tuil ~ ^(u) 8k B TNp 2 L 3 ~ <j> n (u) , . 

X - [UJ) - NMuj ^ ~ vr 4 9 M 3 ca ^ n 4 ' 1 ' 

n=l 71 J n=l 
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One can see that the contributions of higher violin modes are very small due to the factor re -4 in the 
sum. Taking <p n = £ n w and assuming k e L 3> 1, we find the following expression for the violin- mode 
contribution to the thermal-noise spectrum, 



Xv 45 u ' (15) 

Dependence of thermal noise on wire material and suspension 
parameters 



It can be seen from Eqs. (y) and (|15| ) that the thermal noise increases with the area A of the wire 
cross section. Therefore, it is desirable to use wires as thin as possible. However, the wire thickness 
may not be too small since the stress a = T/A in the wire may not exceed the breaking stress a^. In 
fact, the wires are always operated at a fixed fraction of their breaking stress, 

(7 = 0-0 = K(T hr , (16) 

where k is a numerical coefficient. Typical values of k are from 0.3 to 0.5 (it is undesirable to have larger 
values of k because then events of spontaneous stress release will contribute excess noise ||20fl). Thus 
for a given type of the wire material, the cross-section area A should be proportional to the pendulum 
mass M, according to the relation do = Mg/{NA). For a cylindrical wire, one has I = A 2 /2tt. Then 
we obtain 

2 4k B T ( g 3 E V /2 w 



= i=kBTfa? [ * 3 , I (18) 



( 9 3 E 



1/2 



45 rv \2irMN 3 a\ 



L0 



where p v = p/A is the volume mass density of the wire which depends only on the material used. 
All the parameters in Eqs. (fl7|) and ($J$ ) are easily measured except for the wire loss function 



cf> w . A number of experiments were recently performed 14, 15, 21, 22, ^3], |24| to study internal 
losses of various wire materials (e.g., steel, tungsten, fused quartz, and some others). However, the 
exact frequency dependence of the wire loss function (j) w (u>) is not yet completely understood. In 
many experiments (j) w was measured only at few frequencies and experimental uncertainty of results 
was often quite large. Moreover, there are discrepancies between results of different experiments. 
Therefore, it is sometimes difficult to make certain conclusions about the behavior of (f) w (u). 

A well known dissipation mechanism for thin samples in flexure is the so-called thermoelastic 
damping [19|. As a wire bends, one side contracts and heats and the other expands and cools. The 
resulting thermal diffusion leads to the dissipation of energy. The corresponding loss function is 

1 + U> Z T A 

where A is the relaxation strength and f is the relaxation time. The loss function has its maximum 
(ft = A/2 at u = r^ 1 (this is called the Debye peak). This behavior is characteristic for processes in 
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which the relaxation of stress and strain is exponential and occurs via a diffusion mechanism. For the 
thermoelastic damping, one has 

. ETa 2 d 2 , . 

where a is the linear thermal expansion coefficient, C v is the specific heat per unit volume, d is the 
characteristic distance heat must flow, and D is the thermal diffusion coefficient, D = g/C v , where g 
is the thermal conductivity. For a cylindrical wire of diameter d, the frequency of the Debye peak is 

'-5f fi 4 (21) 

For thin metallic wires (d ~ 100 fim) at the room temperature, the Debye peak frequency is typically 
from few hundred Hz to few kHz. Therefore at the frequency range near 100 Hz, we are usually far 
below the Debye peak, and 

4> w (oj) ~ Auf = PAoj, (22) 

where ~ A/(1.3tt 2 L>). 

According to a recent experiment by Huang and Saulson fl5|| , internal losses in stainless steel 
wires are in good agreement with predictions of thermoelastic damping, with 4> w (co) exhibiting the 



characteristic frequency dependence of Eq. (IS). On the other hand, the loss function for tungsten 
wires was nearly constant, increasing slightly at high frequencies (above 500 Hz). <p w for tungsten 
wires increased with the wire cross-section area A, but the exact functional dependence of <j) w on A 
is unclear as only three different wire diameters were examined. In some other experiments, the loss 
functions for various materials were found to be nearly constant over a wide frequency range. In a 
recent experiment by Cagnoli et al. [24|, internal damping of a variety of metallic wires was found to 
be well modelled by the loss function of the form 

(f> w (u) = (f>Q + 0ted(w), (23) 



where </>ted( w ) is the thermoelastic-damping loss function of Eq. (|19| ) and </>o is a frequency-independent 
term. Unfortunately, the dependence of <j)Q on the wire diameter was not examined. It can be assumed 
that the thermoelastic damping is a basic dissipation mechanism, but for some materials it is masked 
by other processes. When those additional losses (whose nature is still a matter of controversy) are 
small, the characteristic frequency dependence of Eq. (|l9|) may be observed. However, when the losses 
due to the thermoelastic damping are very small (which happens, for example, in the case of thin 
Invar and tungsten wires), then additional losses prevail, leading to c/> w which is nearly constant far 
from the Debye peak. 

In what follows we will consider two possibilities: (i) a constant loss function W and (ii) the loss 
function of Eq. (E^) which is characteristic for the thermoelastic damping at frequencies well below the 
Debye peak. We might assume that for some materials the true behavior is somewhere between these 
two extreme variants. For example, for tungsten wires, 4> w is nearly frequency-independent from 50 
to 500 Hz, but still increases to some extent with the wire cross-section area A, as one should expect 
from Eq. (62). 



3.1 A constant loss function 

For a constant <p w , the dependence of the thermal- noise spectrum on various parameters is given 
directly by Eqs. (17) and (|l8|). For the pendulum-mode contribution, we find 
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for constant M and ao, x 2 oc N~ l l 2 ; 

for constant M and N, x 2 oc Cq 1 ; 

for constant A" and ao, x 2 oc M -1 / 2 . 

For the violin-modes contribution, we find 

for constant M and ao, x 2 oc iV -3 / 2 ; 

for constant M and AT, x 2 oc <r^" 3 ; 

for constant N and <to, x 2 oc M -1 / 2 . 

The allowed stress a§ is a property of the wire material (which is also true for E, p v , and 4> w ), so 
changing ao means taking wires made of different materials. Clearly, it is desirable to have a material 
with a large value of do, but what decides is the value of the factor A w = S 1//2 ^ w /o"o for the pendulum 
mode and A w = p^E 1 ! 2 ^ j a\ for the violin modes. The factor A w comprises all the parameters in x 2 
which characterize the wire material. 
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Figure 1: The thermal-noise displacement spectrum y/x 2 (uj) for a multi-loop pendulum suspension 
with tungsten wires: A^ = 4, Q w = 1.3 x 10 3 (solid line); A^ = 16, Q w = 1.3 x 10 3 (dashed line); 
N = 16, Q w = 4.0 x 10 3 (dash-dot line). 

One may see that taking multi-loop suspensions with large numbers of wires may help to reduce 
the thermal noise. As an example, let us consider tungsten wires of the type examined by Huang 



and Saulson [15]. The relevant parameters are E ~ 3.4 x 10 Pa, a^, r ~ 1671 MPa, p v ~ 1.93 x 10 
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kg/m 3 . We also take M = 10.8 kg, L = 45 cm and k = 0.5 (the wires are operated at one half of 
their breaking stress), like in suspensions of the LIGO test masses. According to the data by Huang 
and Saulson fli~5| , the loss function is nearly frequency-independent from 50 to 500 Hz, but depends on 
the wire diameter. For a two-loop suspension (N = 4), the wire diameter should be d ~ 200 /im, and 
the corresponding quality factor Q w = <f>~ 1 can be estimated to be Q w ~ 1.3 x 10 3 . For an eight-loop 
suspension (N = 16), the wire diameter should be d ~ 100 /im, and the corresponding quality factor 
can be estimated to be Q w ~ 4.0 x 10 3 . In Fig. 1 we plot the thermal-noise displacement spectrum 
\fx 2 {u) for the room temperature (T = 295 K) for three possibilities: (a) N = 4, Q w = 1.3 x 10 3 ; 
(b) N = 16, Q w = 1.3 x 10 3 ; (b) N = 16, Q w = 4.0 x 10 3 . We see that for a constant loss function, 
the thermal noise is reduced by increasing the number of wires. The spectral density x 2 {uj) scales as 
N^ 1 / 2 for frequencies near 100 Hz (where the pendulum mode dominates), in accordance with our 
analysis. Also, if the decrease of 4> w with the wire diameter is taken into account, the increase in the 
number of wires is even more helpful. 



3.2 Thermoelastic loss function 

If we take the loss function of Eq. (p?2|), then the thermal-noise spectrum is given by 

2 4k B T ( g 5 EM \ 1/2 1 

/ g 5 EM ' 1/2 



<{ u) = -k B T^ . (25) 



„2(,.\_ ° i_TflJr2 

1 o— \ ■ ■"> ,^ 

The dependence of the thermal-noise spectrum on various parameters can be characterized as follows. 
For the pendulum-mode contribution, we find 



for constant M and <7o, x 2 oc N~ 3 / 2 ; 
for constant M and N, x 2 oc <Jq 2 ; 
for constant and a , x 2 oc 

MV2. 

For the violin-modes contribution, we find 

for constant M and (Tq, x 2 oc N~ 5 ^ 2 ; 

for constant M and N, x 2 oc Gq ; 

for constant A^ and gq, x 2 oc M 1 / 2 . 

Now, the dependence of x 2 on the wire material is given by the factor A w = [3E 1 / 2 ja\ for the pendulum 
mode and A w = fip^E 1 ! 2 j o§ for the violin modes. So, the value of the allowed stress o"o in this situation 
is more important than for the case of constant W . 

One may see that in the case of the thermoelastic damping the thermal noise may be reduced 
to a larger extent by increasing the number of wires, as compared to the case of constant W . As 
an example, let us consider wires made of stainless steel (AISI 302), which were examined by Huang 
and Saulson |0|. The relevant parameters are E ~ 1.9 x 10 11 Pa, dbr — 1342 MPa, p v ~ 8.0 x 10 
kg/m 3 . The losses are dominated by the thermoelastic damping mechanism. Taking a ~ 1.6 x 10 -5 
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1/K, C v ~ 4.8 x 10 6 J/(K m 3 ), g ~ 16.3 J/(K m s) and T = 295 K, one obtains A ~ 3.0 x 10~ 3 
and /? ~ 68.6 s/m 2 . We also take M = 10.8 kg, L = 45 cm and k = 0.5, like in suspensions of the 
LIGO test masses. The thermal-noise displacement spectrum y/x 2 (u>) is plotted in Fig. 2 for three 
possibilities: (a) N = 4 (then d ~ 224 /im and / ~ 176 Hz); (b) N = 8 (then d ~ 159 fim, and / ~ 352 
Hz); (b) N = 16 (then <i ~ 112 /im, and / ~ 703 Hz). The conclusion is that the thermal noise may 
be significantly reduced by increasing the number of wires. The numerical results confirm that the 
proportionalities x 2 oc iV- 3 / 2 and x 2 oc N 5 ' 2 are valid for frequencies well below the Debye peak /. 
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Figure 2: The thermal-noise displacement spectrum \JxHuJ) for a multi-loop pendulum suspension 
with stainless steel wires: N = 4 (solid line); N = 8 (dashed line); N = 16 (dash-dot line). 



3.3 Comparison between different materials 

We would like to compare the thermal-noise performance of a multi-loop suspension for different wire 
materials. For example, the tungsten wires examined by Huang and Saulson |l5j have rather low 
breaking stress of 1671 MPa. There exist tungsten wires with higher breaking stress; for example, 



Dawid and Kawamura [14] experimented with tungsten wires for which they measured o"b r = 2037 
MPa. It would be interesting to compare between tungsten wires with different breaking stress but 
with the same loss function. On the other hand, the comparison between wires made of tungsten and 
stainless steel will clarify how the difference in the loss mechanism (frequency-independent 4> w versus 
the thermoelastic damping) affects the thermal-noise spectrum. To this end, we also would like to 
consider a situation in which wires made of stainless steel have all properties as above except for the 
losses being dominated by a mechanism with frequency- independent </> w , instead of the thermoelastic 
damping. 
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We consider an eight-loop suspension (N = 16), with M = 10.8 kg, L = 45 cm and n = 0.5 (like in 
LIGO), and examine four possibilities: (a) tungsten wires as considered in see Sec 3.1, with o"br = 1671 
MPa (this gives d ~ 100 /im) and Q w = 4.0 x 10 3 ; (b) tungsten wires with different breaking stress, 
dbr = 2037 MPa (this gives d ~ 91 jiva) and the same quality factor, Q w = 4.0 x 10 3 ; (c) stainless 
steel wires as considered in Sec 3.2 (<7b r = 1342 MPa, d ~ 112 /im), with the thermoelastic damping 
mechanism (/3 ~ 68.6 s/m 2 , / ~ 703 Hz); (d) stainless steel wires with the same parameters, but 
with a frequency-independent loss function, Q = 2.0 x 10 3 (this value is close to the one given by the 
thermoelastic damping near 120 Hz). The resulting thermal- noise displacement spectra ^x 2 (u) are 
shown in Fig. 3. One can see that the violin resonances of the stainless steel wires appear at higher 
frequencies (due to smaller density). On the other hand, the tungsten wires exhibit smaller thermal 
fluctuations at the frequency range between 50 and 200 Hz. The thermal noise is reduced by using 
wires with larger breaking stress, provided the other parameters remain the same. 
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Figure 3: The thermal-noise displacement spectrum yJx 2 {oj) for an eight-loop pendulum suspension 
(N = 16): tungsten wires with = 1671 MPa and Q w = 4.0 x 10 3 (solid line); tungsten wires with 
<7br = 2037 MPa and Q w = 4.0 x 10 3 (dashed line); stainless steel wires with <7br = 1342 MPa and 
thermoelastic damping (dash-dot line); stainless steel wires with a hT = 1342 MPa and Q = 2.0 x 10 3 
(dotted line). 

3.4 Optimization of the pendulum length 

The thermal-noise spectrum depends on the pendulum length L. For frequencies well below the first 
violin resonance, lo 2 CwJ, the pendulum- mode contribution dominates and the spectral density x 2 (to) 
is proportional to L~ 2 . However, by increasing L, one not only decreases the thermal fluctuations 
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due to the pendulum mode, but also brings the violin resonances to lower frequencies, as u n oc L . 
This effect is illustrated in Fig. 4, where the displacement spectrum yJx 2 {uS) is shown for an eight- 
loop suspension with stainless steel wires of various length. (We take M = 10.8 kg, k = 0.5, and 
stainless steel wires with properties listed in Sec. 3.2.) Due to this competition between two opposite 
tendencies, the choice of the pendulum length is a delicate matter which depends on where in the 
spectrum the seismic perturbations and the photon shot noise prevail over the thermal fluctuations 
and on properties of expected gravitational-wave signals. 
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Figure 4: The thermal-noise displacement spectrum ^x 2 (u>) for an eight-loop pendulum suspension 
(JV = 16) with stainless steel wires of various length: L = 30 cm (solid line); L = 45 cm (dashed line); 
L = 60 cm (dash-dot line). 



4 Discussion 

Our analysis brings to an observation that the thermal noise in pendulum suspensions can be sig- 
nificantly reduced by using multi-loop configurations with a large number of wires. However, before 
implementing this conclusion one should consider a number of issues. First, our analysis is valid only 
if the losses are dominated by the internal friction in the pendulum wires and all other sources of 
dissipation are made negligible by careful experimental design. However, as was shown recently by 
Huang and Saulson [p^] , the sliding friction in the suspension clamps is often important as well. If 
this is the case, a large number of suspension loops will only sever the dissipation and thereby increase 
the thermal fluctuations. Therefore, if one wants to use multi-loop suspensions, a special care should 
be paid to the design of clamps. Another technical problem is to make a suspension in which all the 
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loops will be equally loaded. One more issue which should be carefully studied is the effect which may 
have a large number of suspension wires on the internal resonances of the suspended test mass. 
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